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The influence of the nuclear medium on the production of charged hadrons in semi- 
inclusive deep-inelastic scattering has been studied by the HERMES experiment at DESY 
using a 27.5 GeV positron beam. The differential multiplicity of charged hadrons and 
identified charged pions from nitrogen relative to that from deuterium has been measured 
as a function of the virtual photon energy v and the fraction z of this energy transferred 
to the hadron. There are observed substantial reductions of the multiplicity ratio R\ { at 
low v and at high z, both of which are well described by a gluon-bremsstrahlung model 
of hadronization. A significant difference of the ^-dependence of R^ is found between 
positive and negative hadrons. This is interpreted in terms of a difference between the 
formation times of protons and pions, using a phenomenological model to describe the v- 
and z-dependence of R\ { . 



PACS numbers: 13.60.Hb, 13.60.Le, 13.60.-r, 24.85.+p 
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In deep-inelastic scattering (DIS) an incident 
lepton interacts with a target via the exchange 
of a gauge boson between the lepton and a par- 
ton (a quark or an antiquark). The struck par- 
ton is subjected to a sequence of hard parton- 
and soft hadron-production processes, resulting 
in the formation of hadrons. By carrying out 
DIS experiments on nuclear targets, it is possi- 
ble to study the hadronization process during the 
time period immediately after the quark has been 
struck by the virtual photon pj|. In the simplest 
scenario the nucleus, which has the size of a few 
fm, acts as an ensemble of targets with which the 
struck quark or the produced hadron may interact. 
If an interaction occurs, the number of leading 
hadrons produced per DIS event and per nucleon 
is reduced compared to that for a free nucleon. 
The reduction of hadron multiplicity depends on 
the distance traversed by the struck quark be- 
fore the hadron is formed, the (unknown) quark- 
nucleon cross section and the (known) hadron- 
nucleon cross section. Hence, measurements of 
the multiplicity of hadrons produced on nuclei can 
provide information on the space-time structure of 
the hadronization process. 

In this paper we present the results of deep- 
inelastic positron scattering measurements on 14 N 
and 2 H targets carried out by the HERMES Col- 
laboration at DESY. The reduction of the hadron 
multiplicity on 14 N relative to that on 2 H has been 
determined from the data. Deuterium was used as 
a reference target instead of hydrogen to account 
for the difference in scattering from neutrons or 
protons. 

In the past, semi- inclusive leptoproduction of 
hadrons from nuclei was studied at SLAC with 
electrons f| and at CERN and FNAL with high- 
energy muons by EMC || and E665 (|. In these 
experiments the effect of the nucleus on the early 
stages of hadronization appeared in the depen- 
dence of the multiplicity ratios on the energy v 
of the virtual photon (as defined in the target rest 
frame). The effect is most prominent if v ranges 
from a few GeV to a few tens of GeV. The z^-range 
covered by the present experiment (7 - 23 GeV 
was used in the present analysis) is low enough to 
be optimal for studying the influence of the nu- 
clear medium on the hadronization process, while 



it is still high enough to be in the scaling regime, 
i.e. for the DIS-framework to apply. 

The theoretical description of hadronization is 
difficult, as it is impossible to treat the problem 
in a perturbative QCD-framework. Although the 
nuclear environment can be used to study certain 
aspects of the hadron formation process, the same 
difficulties apply to the description of nuclear ef- 
fects on hadronization. Therefore, one has to re- 
sort to model calculations. The HERMES data 
on hadron production from nuclei are compared 
to the gluon-bremsstrahlung-model calculations of 
Ref. jjj, and the phenomenological intranuclear 
reinteraction model ofRefs. ||[7j. 

In the gluon-bremsstrahlung model, the struck 
quark is assumed to lose energy via the emission 
of gluons until, in the case of meson formation, a 
qq configuration is formed consisting of the struck 
quark and an antiquark originating from the last 
emitted gluon |||| [rJ- If this last step occurs in- 
side the nucleus, the nuclear environment affects 
the hadron multiplicity because the meson inter- 
acts with the nuclear medium with a sizable cross 
section. Moreover, the interaction of the initial 
quark with the nuclear medium causes the emis- 
sion of additional soft gluons. On the other hand, 
the initial — possibly small — qq configuration 
represents a color dipole which may have a re- 
duced probability of interaction with its environ- 
ment (Color Transparency ]ll|,0]). An estimate 
of the combined effect of the soft gluon radiation 
and Color Transparency on the multiplicity ra- 
tio in the framework of the gluon-bremsstrahlung 
model [|j yields 2-3% for the kinematics of the 
present experiment. 

An important parameter characterizing the 
hadronization process is the formation time (in 
the target rest frame) t'j, which is defined as the 
mean time elapsed between the moment the quark 
is struck and the creation of the leading hadron 
h. In the gluon-bremsstrahlung model ||, the 
probability distribution of t'j may be evaluated 
for given values of z and is, where z represents 
the fraction of the energy v transferred to the 
hadron. The mean value of tj exhibits an approx- 
imate (1 — ^independence. In the phenomenolog- 
ical models ||[7j , on the other hand, the formation 
time is assumed to be equal to the product of a 
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time constant intrinsic to the kind of hadron 
produced, and a time-dilation factor v jm where m 
represents a mass. In the past, several functional 
forms for ih have been considered j6|,^| JTl|Jl^ |l5| 
for instance v) — t^zv /rrih or vj — Thv/m q , where 
mh (jnq) represents the mass of the hadron (con- 
stituent quark). Alternatively, one may take for 
m the kinematic expression derived by Berger fl6|| 
for the mass of the offshell struck (current) quark. 
Under the assumption that the mass of the initial 
quark and that of the produced hadron can be ne- 
glected, one finds th — Th\J z{l — z)/p'j,u, where 
Pt represents the component of the hadron mo- 
mentum transverse to the direction of the virtual- 
photon momentum. 

The experimental results are presented in terms 
of the multiplicity ratio R\ t , which represents the 
ratio of the number of hadrons of type h produced 
per DIS event for a nuclear target of mass A to 
that from a deuterium target (D): 



R h M {z,v) 



N£{z,v) I N»{z,v) 



N e{ v ) 



N?{v) 



(1) 



with Nhiz^) the number of semi-inclusive 
hadrons in a given (z, ^)-bin, and N e (y) the num- 
ber of inclusive DIS positrons in the same v-h\n. 
For the purpose of the present analysis, the mul- 
tiplicity ratio was determined as a function of v 
and z, while integrating over all other kinematic 
variables. The data for R 1 ^ are only weakly de- 
pendent on either Q 2 or p\ ]2j| . 

The HERMES experiment at DESY makes use 
of the 27.5 GeV positron storage ring of HERA. 
The 14 N and 2 H target gases are injected into 
a tubular open-ended storage cell through which 
the beam passes. The cell provides a 40 cm long 
target with areal densities of up to 6 x 10 15 
nucleons/cm 2 for 14 N. The dead time of the data 
acquisition system was observed to be less than 
5% even at the highest luminosities of about 10 33 

— 2 —1 

cm s . 

In the HERMES spectrometer both the 
scattered positrons and the produced hadrons are 
detected and identified within an angular accep- 
tance of ± 170 mrad horizontally, and ± (40 - 
140) mrad vertically. The trigger for scattered 
positrons was formed from a coincidence between 
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FIG. 1. Scatter plot of the hadron (or pion) energy 
Eh {Ett) and the energy transfer v. Lines representing 
constant values of z are shown as well. 

several scintillator hodoscope planes and a lead- 
glass calorimeter. The trigger required an en- 
ergy of more than 3.5 GeV deposited locally in 
the calorimeter. Positron identification was ac- 
complished using the calorimeter, the second ho- 
doscope, which functioned as a preshower counter, 
a transition-radiation detector, and a threshold 
gas Cerenkov counter filled with a mixture of N2 
and C4-F10 at atmospheric pressure. This system 
provided positron identification with an average 
efficiency of 99 % and a hadron contamination of 
less than 1 %. Hadrons with an energy Eh > 1.4 
GeV could be distinguished from leptons. Fur- 
thermore, in the momentum range between 4 and 
13.5 GeV, pions were identified with the help of 
the Cerenkov counter. 

Events were selected by imposing constraints on 
the four- momentum squared — Q 2 of the virtual 
photon, the invariant mass of the photon-nucleon 
system W, the Bjorken scaling variable x — 2 ^ - 
(with m p the proton mass), and the energy frac- 
tion of the virtual photon y = v / ' Eh ea m- For each 
event it was required that x > 0.06, Q 2 > 1 GeV 2 , 
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W > 2 GeV, and y < 0.85. The requirements on 
W and y are applied to exclude nucleon resonances 
and to limit the magnitude of the radiative correc- 
tions, respectively. The constraints on (x, Q 2 ) are 
applied to exclude the kinematic region where an 
anomalous ratio of the longitudinal to transverse 
cross sections for inclusive deep-inelastic scatter- 
ing from 14 N has been observed |D| . The anomaly 
can be interpreted as being due to the absorp- 
tion of the virtual photon by a correlated quark 
pair or by a nuclear meson pOfl . Either inter- 
pretation breaches the assumption of incoherent 
lepton-quark scattering, inherent in deep-inelastic 
scattering. 

The instrumental threshold on Eh is larger for 
pions than for hadrons. This implies, as is illus- 
trated in Fig. [l], that e.g. the pion z-acceptance is 
restricted to rather large values as v decreases. 
Hence, to ensure that the ^-dependence of the 
multiplicity ratio does not correspond to a strong 
variation of the mean z-value, the presented data 
are confined to z > 0.2 and v > 7 GeV for hadrons 
and z > 0.5 and v > 8 GeV for pions. 

Under the kinematic constraints described 
above, the number of selected DIS events is 0.88 
(1.05) xlO 6 for 14 N ( 2 H) if v > 7 GeV. For v > 8 
GeV these numbers decrease to 0.76 xlO 6 for 14 N 
and 0.91 x 10 6 for 2 H. The number of hadrons with 
z > 0.2 equals 227 (288) xlO 3 , and the number of 
pions with z > 0.5 is 36 (48) xlO 3 for 14 N ( 2 H). 

The HERMES data have been corrected for ra- 
diative processes using the codes of Refs. pl| , p2| . 
As is commonly done |^|,^| , the contributions from 
nuclear elastic scattering, quasi-elastic scattering 
and DIS were treated independently. Coherence- 
length effects 23 1 are not included in this ap- 
proach. The code of Ref. [@ was modified to 
include the measured semi-inclusive DIS cross sec- 
tions for 14 N. The size of the radiative correction 
was found to be negligible in most of the kinematic 
range, with a maximum of 3% at the highest value 
of v. The correction is so small because most of 
the contributions to the radiative corrections can- 
cel in the multiplicity ratio. 

The systematic uncertainty of the present data 
is 3% or less. The uncertainty is small due to the 
fact that double ratios of semi-inclusive and inclu- 
sive yields are measured. The main contributions 
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FIG. 2. Charged hadron multiplicity ratio R%[ as 
a function of v for values of z larger than 0.2. In 
the upper panel the CERN |] and SLAC [§ data for 
Cu are compared to various phenomenological calcu- 
lations taken from the original publications In 
the lower panel the HERMES data on 14 N are repre- 
sented by solid squares, while the open star represents 
the CERN data point on 12 C and the open square 
the SLAC data point on 12 C. The error bars represent 
the statistical uncertainty only. The systematic un- 
certainty of the HERMES data is < 3%. The curves 
are described in the text. 

to the systematic uncertainty arise from radiative 
corrections (< 2%), overall efficiency (< 1.5%), 
and diffractive p°-meson production (<1.5%) [P5| . 
It has been verified that the acceptance for semi- 
inclusive hadron production is the same for both 
targets, by plotting the multiplicity ratio versus 
hadron angle and y\. These ratios are constant 
to within the expected statistical precision except 
for the data beyond p|i = 1 GeV 2 , where an en- 
hancement similar to that reported by Ref. || 
is observed. Because of the steep decline of 
semi- inclusive hadron production with p 2 ^, this 
enhancement does not contribute to the present 
data. The nuclear attenuation of pions resulting 
from the decay of p-mesons formed in the frag- 
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mentation process (which constitute about 20% 
of the semi- inclusive pion yield) was found to 
have a (/-dependence similar to that of the direct- 
fragmentation pions. No correction for this pro- 
cess was therefore applied. 

The results for the multiplicity ratio for all 
charged hadrons with z > 0.2 are presented as 
a function of v in Fig. |^ together with data of pre- 
vious experiments. In the top panel, existing data 
on copper are shown together with the origi- 
nal phenomenological calculations that were used 
to interpret the data. In the lower panel of Fig. ^, 
the present data on 14 N are displayed together 
with data of previous experiments at CERN || 
and SLAC || on 12 C. The SLAC data point is 
lower than the present data by five standard devi- 
ations. For the SLAC data no systematic uncer- 
tainty is available from the original publication. 
Moreover, no corrections were made for the target- 
mass dependence of the inclusive DIS cross sec- 
tion, an effect that was not known at the time of 
the analysis. Such a correction is needed because 
a semi-inclusive cross-section ratio was measured 
at SLAC instead of the multiplicity ratio defined 
in Eq. (1). An estimate of the correction for this 
effect - based on the HERMES 14 N data - results 
in a 4% increase of i?5u, which reduces the dis- 
crepancy by a factor of two. 

The present data for R^j are observed to in- 
crease with increasing v and approach unity for 
high v, which is consistent with the EMC data. 
This behaviour suggests that for large values of 
hadron formation takes place mainly outside the 

TABLE I. Parameters used in the phenomenologi- 
cal one and two time-scale model calculations shown 
in Figs. 2 and 3. In the left two columns the line 
type used in the figures and the expression used for 
the formation time are listed, while the three columns 
on the right give the values of the three cross sections 
(in mb) used in the calculations. 
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FIG. 3. The multiplicity ratio as a function 
of z for all charged pions (open circles) and all 
charged hadrons including pions (closed squares) . The 
full curve represents a gluon-bremsstrahlung model 
calculation for pions. The dotted, dashed and 
dot-dashed curves represent phenomenological forma- 
tion-time calculations. 

14 N nucleus, and that the interaction of the struck 
quark with the nuclear medium is weak. 

For the interpretation of the data it is useful 
to provide some details on the calculations from 
Ref. (D shown in the upper panel of Fig. ||. Similar 
calculations are also used when discussing the z- 
dependence of the 14 N data below. The long-dash 
dotted curve represents a calculation in the frame- 
work of the one time-scale model Q , in which it is 
assumed that the influence of the nuclear medium 
on the hadron multiplicity is due to quark- nucleon 
scattering for t < t h f and hadron-nucleon scatter- 
ing for t > th, with cross sections a q and ah re- 
spectively. In the calculation shown t h j = Chzv 
with Ch — 1 fm/(GeVc) has been assumed, 
and the parameters listed in table | have been 
used. The other curves in fig. |^ are three ex- 
amples of two time-scale model calculations, in 
which the formation time is assumed to be given 
by tf = (l — \n(z))zv / (hic) with c the speed of 
light and k ~ 1 GeV/fm the string tension fTUnjl- 
In these models an additional parameter a s is in- 
troduced that represents the interaction between 
the open string and the medium occurring after a 
time interval r c = tj ~ zv/(kc). The expressions 
for and r c have been derived from the color- 
string model Q. The parameters used for the 
three calculations shown are also listed in table 0. 



G 



The long-dashed curves give a fairly good descrip- 
tion of both the Cu and the present 14 N data. 
However, it should be noted that a two time-scale 
model calculation with a s = is equivalent to 
a one time-scale model calculation except with a 
different expression for v\. For this reason we will 
use only the one time-scale model when present- 
ing the final parameterization of the data near the 
end of this paper. 

In order to further investigate the various possi- 
ble descriptions of the hadron formation process, 
the z-dependence of was extracted from the 
present measurements for both hadrons and iden- 
tified pions. The multiplicity ratio for hadrons 
(pions) with v > 7 GeV (8 GeV) is presented as a 
function of z in Fig. |[ 

No significant difference between data for 
hadrons and pions is observed here. The de- 
crease of R'lj with z, as observed at large z for 
the first time by the present experiment, is at 
variance with the phenomcnological-model calcu- 
lations, which predict an increase of R h M with z. 
This is demonstrated by the long-dash-dotted (one 
time-scale) and long-dashed (two time-scale) cal- 
culations shown in Fig. ^. The curves have been 
obtained in the same way as the corresponding 
curves in Fig. ^ only using the 14 N matter density 
instead of the one for Cu. It is concluded that 
the (z, ^-dependence of the formation time is not 
given by the simple expressions mentioned 
above. It is significant that the expression for 
implied by the kinematic z-dependence of the in- 
variant mass of the current quark also fails 
to reproduce the observed dependence of R^ on 
either z or p| 

The solid curve, on the other hand, gives a 
good account of the data. It represents the re- 
sult of a calculation for RJ f within the gluon- 
bremsstrahlung model M. Note that the calcu- 
lation applies to pions only. The basic mechanism 
that causes the decrease of R^ at large values of z 
is the following. A hadron with large z originates 
from a quark that has emitted only few gluons 
of relatively small energy; otherwise its z-value 
would have been lower. As the bare quark contin- 
uously emits gluons, the emission of only few glu- 
ons corresponds to a small formation time. Hence, 
the probability that the high-z hadron is subject 
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FIG. 4. The multiplicity ratio as a function of v 
for charged pions with z > 0.5. The solid curve repre- 
sents a gluon-bremsstrahlung model calculation. The 
dotted curve is the result of a one time-scale model 
calculation assuming a (1 — z)i/-dependence of the for- 
mation time. 

to rescattering is largest, explaining the decrease 
of R^ at high z. 

In order to be able to study phenomenological 
concepts such as the quark-nucleon cross section 
dq and the formation time t 1 }, the z-dependence 
of R\,f has been parameterized in the framework 
of the one time-scale model. Using the expres- 
sion tj — c/j(I — z)v ■, as suggested by the gluon- 
bremsstrahlung model M, the dotted curve in 
Fig. H has been obtained. It represents the result 
of a fit with the factor Ch a free parameter, and 
Oq and ah set to the values listed in table [jj. It is 
concluded that a phenomenological description of 
the z-dependence of the R\ t data can be obtained 
if an ad-hoc (1 — z)-dependence of the formation 
time is assumed. 

The //-dependence of the multiplicity ratio for 
charged pions with z > 0.5 is shown in Fig. ||. 
The calculation with the gluon-bremsstrahlung 
model H is in fair agreement with the data. A 
similar description of the data is obtained using 
the one time-scale model with tj — c n (l — z)v, as- 
suming a q — mb and <ih = 25 mb (dotted curve 
in Fig. ^J). In a two-parameter fit with and 
<7 q treated as free parameters, the resulting value 
for a q is consistent with zero with an uncertainty 
of only ±0.4 mb (or ±0.2 mb if the hadrons are 
considered) . It is concluded that in the context of 
the phenomenological models, the interaction of 
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the struck quark with the nuclear medium is in- 
deed very small, which is in accordance with the 
results of Refs. [s|Jl5|]. Hereafter, we extract fur- 
ther information about the hadron formation time 
within the context of the one time-scale model 
with o q rs 0. 

The results presented thusfar concern the sum 
of positive and negative hadrons. For both pi- 
ons and hadrons the charge states can be studied 
separately. In the upper (lower) panel of Fig. |j| 
the multiplicity ratios for positive and negative 
hadrons (pions) are displayed as a function of v. 
The same requirements on the variable v as before 
have been imposed: v > 7 GeV for the hadron 
data and v > 8 GeV for the pion data. However, 
for the purpose of this comparison and in order to 
increase the pion statistics, the requirement z > 
0.2 has been imposed in both cases. 

The data show that the multiplicity ratio R\j 
is the same for positive and negative pions. This 
observation — together with the fact that the 
tt^ 1 — N interaction cross sections are about equal 
- is consistent with the likely assumption that 
the formation times for the two charge states are 
the same. 

In contrast, a significant difference is observed 
between R t ^ 1 of positive and negative hadrons. 
Since the multiplicity ratios of positive and nega- 
tive pions are measured to be equal, this difference 
must be caused by other hadrons, such as protons, 
antiprotons and kaons. A Monte Carlo study with 
the LUND model |2(| has been performed in order 
to estimate the relative yield of various hadrons 
in the acceptance of the spectrometer. The study 
shows that the contribution of antiprotons (pro- 
tons) to the negative (positive) hadron sample, 
averaged in the z/-range of 3 - 23 GeV, equals 6% 
(24%). In the same kinematic range, the negative 
and positive kaon contributions are 10 and 13%, 
respectively. Since the effects of the interaction 
between the struck quark and the nuclear medium 
are found to be small, only two effects may ac- 
count for the observed difference between Rzf and 
R 



M 



the difference in the cross sections |27]] for 
rescattering of antiprotons (~ 60 mb), protons (« 
40 mb), negative kaons (ss 23 mb) and positive 
kaons (w 17 mb), or possible differences between 
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FIG. 5. Multiplicity ratios for hadrons includ- 
ing pions (top panel) and identified pions (bottom 
panel) as a function of v. The open (closed) squares 
represent the positive (negative) hadrons. Identified 
pions are represented by open (positive) and closed 
(negative) circles. The curves are parameterizations 
of the data using the one time-scale model assuming 
t) =c h (l-z)v. 

the formation times of baryons and mesons. Not- 
ing that the effective hadron-nucleon cross sec- 
tions, obtained as the weighted sum of individual 
hadron cross sections, are practically the same for 
the positive and negative hadron samples, the ob- 
served difference of R^ and R 1 ^ is most likely 
due to differences in the formation times of pi- 
ons, kaons and/or (anti)protons. Since the proton 
fraction in the positive hadron sample is four times 
larger than the antiproton fraction in the negative 
hadron sample, it is concluded that protons have a 
significantly larger formation time than pions. In 
order to substantiate this conclusion, additional 
data with good pion, kaon and proton identifica- 
tion are needed. 

In the absence of a theoretical model predict- 
ing the observed difference between R\ d and R\j , 
the data of Fig. ^ have been parameterized using 
the one time-scale model 101. In these calcula- 
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tions, the formation time is assumed to be given 
by tj = c/j(l — z)v to ensure a proper descrip- 
tion of the z-dependence of R^ , and an effective 
hadron-nucleon cross section of 25 mb is used for 
all hadrons. Initially both a q and Ch were used 
as free parameters. In all four cases (h + , h~, 
7r + , and tt~ ) the fits gave a good account of the 
data, while the fitted value of a q was found to be 
consistent with zero within uncertainties ranging 
from 0.3 mb for 7r + to 1.2 mb for h _ . This im- 
plies that the dominant effect of the nucleus on 
R\j is due to hadron-nucleon interactions. Set- 
ting Gq equal to zero, the fits were repeated yield- 
ing once more a good description of the data (see 
Fig. ||) and the following set of fit parameters: 
c w - = 1.38±0.21 fm/(GeVc), = 1.37±0.18 
fm/(GeVc), c h - = 1.32±0.16 fm/(GeVc), and c h + 
= 3.49±0.51 fm/(GeVc). The quoted uncertain- 
ties equal the quadratic sum of the statistical er- 
rors of the fit parameters and the variations of 
these parameters when the hadron-nucleon cross 
section is varied by ±2 mb. The fitted values of 
Ch quantitatively confirm our earlier conclusions 
concerning the equality of the formation times of 
positive and negative pions, and the large differ- 
ence between the formation times of positive and 
negative hadrons. 

In summary, high accuracy data on the at- 
tenuation of charged hadrons and identified pi- 
ons in deep-inelastic scattering on 14 N relative to 
2 H have been obtained by the HERMES experi- 
ment. The ^-dependence of the multiplicity ratio 
of hadrons (h + and h _ ) has a similar behaviour 
as the CERN muon H data. In the previously 
unexplored high z-region covered by the present 
experiment, a decrease of R\ 4 with z is observed. 
This observation is in agreement with the predic- 
tion of the gluon bremsstrahlung model ||, but 
is at variance with the prescription often used in 
phenomenological models that the hadron forma- 
tion time is proportional to z, and with the z- 
dependence implied by Ref. |ll|. In contrast to 
what is observed for positive and negative pions, 
a difference is found in the ^-dependence of R^ 
for positive and negative hadrons. This difference 
is interpreted to indicate that a proton has a larger 
formation time than a pion. 

Additional measurements of hadronization in 



various heavy nuclei with pion, kaon and pro- 
ton identification are needed to clarify the issues 
raised by the present data concerning formation- 
time differences of different kinds of hadrons. Such 
measurements are underway at HERMES. 
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